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RESPONSE OF A RADIAL-BLADED CENTRIFUGAL PUMP TO SINUSOIDAL 
DISTURBANCES FOR NONCAVITATING FLOW 
by Douglas A. Anderson, Robert J. Blade, and William Stevans 

Lewis Research Center 

SUMMARY 

A radial -bladed centrifugal pump was run in water with sinusoidal fluctuations of 
pressure and flow rate imposed at the pump inlet. The pump response was determined 
by measurements of appropriate fluctuating pressures and flow rates. Since the flow 
was noncavitating, zero gain was assumed in computing pump impedance. Furthermore, 
it was shown that the gain computed from the data was zero within the accuracy of the 
measurements . 

The pump impedance was found to be a complex value composed of resistance and 
inertive reactance. The reactance was greater than the resistance except at relatively 
low frequencies (below 3 to 9 Hz, depending on the flow rate). The resistances com- 
puted from the steady -state performance curves agreed reasonably well with the dynamic 
values at zero frequency. Both resistance and inertance were frequency dependent, with 
the greatest changes occurring below 20 hertz. For the impeller -volute combination, 
the resistance increased with frequency while the inertance decreased. For the outlet 
diffuser, the resistance and inertance were almost independent of frequency near design 
flow, but had significant low-frequency trends at higher flows. The high-frequency dif- 
fuser inertance agreed with the value computed from the diffuser geometry. 

An electric circuit model was developed in order to explain the trends of impedance 
with frequency. The model does not require any dynamic measurements . The model 
results agreed reasonably well with the measured data near design flow and at high fre- 
quencies for higher than design flows. However, the model does not explain some of the 
low-frequency -data trends at the higher flows. These trends are probably due to the ef- 
fects of the uneven flow distribution in the impeller and the volute. 


INTRODUCTION 


In launch vehicles, the longitudinal structural oscillations often referred to as 



"pogo” are a serious problem. These oscillations result from the coupling of the ve- 
hicle structure and the vehicle propulsion system in a closed -loop feedback arrangement. 
To analyze these oscillations requires that each component of the closed loop be repre- 
sented by an appropriate transfer function. One of the problems encountered in pogo 
studies is in predicting with sufficient accuracy the response of the propellant pumps to 
variations in inlet pressure and flow rate. 

The pump is commonly represented in terms of gain and impedance as is explained 
in reference 1. While gain and impedance can be related to the steady -state operating 
characteristics of the pump, reference 1 indicates that it is preferable to base the rep- 
resentation on dynamic data. Two attempts to determine pump dynamic response exper- 
imentally are reported in references 2 and 3. Reference 2 presents some dynamic mea- 
surements made in liquid oxygen and liquid hydrogen, and suggests a different repre- 
sentation of the pump than reference 1. Reference 3 gives some results in which the in- 
put perturbations were created by shaking the pump. 

In the investigations to date, the analysis of cavitating dynamic data requires some 
prior assumptions to be made about pump impedance or gain. An experimental investi- 
gation was carried out at the NASA Lewis Research Center in order to obtain more com- 
plete dynamic data. One of the purposes of this investigation was to allow the determi- 
nation of unambiguous values of pump gain and impedance. Reference 4 describes the 
test facility used in this investigation and the design philosophy on which it is based. 

Some work was also done to define the trends of noncavitating pump impedance with 
frequency. This is the main concern of the present report. Data are presented which 
were taken at noncavitating inlet pressures over a range of frequencies and flow rates . 
For noncavitating flow, the dynamic equations become less complicated since gain is 
zero. This allows the pump impedance to be calculated in a more straightforward man- 
ner, and makes the results less sensitive to measurement errors. 

TEST FACILITY AND INSTRUMENTATION 

The test facility used in this investigation is a closed-loop water tunnel. Only a 
brief description of the test facility and instrumentation will be given herein, as a de- 
tailed description appears in reference 4. The test facility has subsystems which are 
used to generate and measure sinusoidal perturbations of pressure and flow rate. A 
schematic diagram of the test facility is shown in figure 1(a), and a photograph in fig- 
ure 2. 
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Test Pump 


The pump selected for these tests has an unshrouded radial -bladed centrifugal im- 
peller and a single discharge volute. A sketch of the pump internal geometry is shown 
in figure 1(b), and a photograph of the impeller in figure 3. The complete steady -state 
performance of the pump impeller is given in reference 5. 

Flow and Pressure Perturbing System 

The purpose of this system is to generate small sinusoidal flow and pressure pulsa- 
tions at the pump inlet. The system is located between the inlet isolation tank and the 
inlet line. One component of the system is a perturbation throttle designed especially 
for the test facility . A precisely controlled variation in throttle flow area produces the 
desired flow and pressure fluctuations. 


Dynamic Pressure Measurements 

These measurements were made with flush -diaphragm, strain-gage pressure sen- 
sors. These sensors are placed at 10 locations, labeled 1 to 10 in figure 1(a). The 
electrical signals from the sensors were measured in phasor form by using a transfer 
function analyzer. Phase data were obtained by comparing the signal to be measured to 
a reference signal (the same signal as was used to drive the perturbation throttle). 


Inlet Dynamic Flow Measurement 

The perturbed flow approaching the pump was computed from the pressure measure- 
ments on the inlet line (stations 1 to 5 of fig. 1(a)). The perturbed flow rate is related to 
these pressures by the acoustical wave equation as given in reference 4 (eq. (1)). 


Vibration Measurement 

The longitudinal motion of the inlet line near the pump was measured with two dif- 
ferent motion sensors. This motion is taken into account in the calculation of the per- 
turbed flow into the pump. Reference 4 showed that if the vibration correction were ne- 
glected, gross errors would result. 
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Outlet Dynamic Flow Measurement 


The perturbed flow leaving the pump was determined by measuring the perturbed 
pressure drop across a specially designed multihole orifice plate. Figure 1(a) shows 
the location of the orifice plate and of the pressure sensors (labeled 9 and 10). The ori- 
fice plate is shown in figure 4. The purpose of the special design is to minimize the 
inertive component of the orifice plate dynamic impedance, as discussed in reference 4. 


Open Orifice and Blocked Orifice 

One of the requirements of the test program was that data be taken at two different 
system exit impedances. The two impedances were obtained by using two different out- 
let orifice configurations, referred to as the open orifice and the blocked orifice. Open 
orifice denotes the multihole orifice plate as originally fabricated. Blocked orifice de- 
notes the same orifice plate with 20 percent of the holes filled with epoxy plastic. The 
blocked orifice had a dynamic resistance 63 percent higher than the open orifice. 


Steady-State Measurements 

Steady-state pressure measurements were made with strain-gage pressure sensors. 
Measurements were made of pump inlet and outlet total pressure, and of wall static 
pressures at stations 6, 8, 9, and 10 of figure 1(a). Steady -state flow rate was mea- 
sured by using a venturi and a bi-fluid manometer . The rotational speed of the pump 
was measured by using a magnetic pickup. 


APPROACH 

Pump dynamic response is generally described in terms of gain, impedance, and 
compliance. This report is concerned primarily with the experimental determination of 
pump impedance. The following discussion first presents the method developed in ref- 
erence 4 for determining impedance and gain. After that, the simplifications possible 
for the case of noncavitating flow are discussed. These simplifications are of special 
interest in this report, which analyzes only noncavitating data. 

The standard equation used to represent pump dynamic response is (ref. 1) 

P D = (m + l)P g - ZQ out (1) 


4 



Equation (1) is a phasor equation in which Pp refers to the perturbed discharge total 
pressure, Pg the perturbed suction total pressure, and Q out the perturbed outlet flow 
rate. (All symbols are defined in appendix A. Throughout this report, pressures and 
flow rates are perturbed parameters unless otherwise stated.) This equation is based 
on the assumption that the perturbed flow is equal to Q out everywhere within the pump. 
Any compliance due to cavitation or to pump structural flexibility is assumed to act at 
the pump inlet but not within the pump itself. In using equation (1), the values of P^, 
Pg, and Q out are obtained from the test data of this investigation. Pump gain m and 
pump impedance Z are then to be calculated from these values . 

In general, the values of m and Z cannot be determined from a single set of Pp, 
Pg, and Q out values. There are two unknowns but only one equation. In the method of 
reference 4, two sets of data are taken at the same steady-state conditions but with dif- 
ferent exit impedances in the test loop. Changing the exit impedance changes the ratios 
of Pj-)/Q ou t and Pg/Q Qut without changing the pump gain m or the pump impedance Z 
Thus, the two sets of data provide two equations from which the two unknowns m and Z 
can be determined. In the test facility, two different exit impedances were obtained by 
using two different configurations of the multihole orifice plate. 

The gain m accounts for the effect of cavitation on the pump response to pressure 
perturbations. For steady -state operation, gain is equal to the slope of the pump head 
rise plotted against inlet head. For dynamic operation, gain will, in general, be a com- 
plex number. For noncavitating flow, gain should theoretically be zero. All the data 
shown in this report are noncavitating, and the zero-gain assumption was used in calcu- 
lating impedance. The significance of the gain data is examined in the section RESULTS 
AND DISCUSSION and in appendix C. It was confirmed that gain is zero within the mea- 
surement accuracy of the data. 

Assuming that the gain is zero, a value of pump impedance can be computed from a 
single set of data. Denoting this value as % nc , 


Z 


nc 



( 2 ) 


The impedance will generally be a complex number. 

Two sets of noncavitating data were taken using two different multihole orifice plate 
configurations. The two sets are referred to as the open-orifice data and the blocked- 
orifice data. The pump impedance can be computed from both sets of data by using 
equation (2); comparing the results then gives an indication of data consistency. The 
pump impedance can also be computed by a simultaneous solution using equation (1) . 
However, for the test data obtained in this investigation, the effect of measurement 
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errors was magnified when the simultaneous solution was used. Consequently, the most 
accurate values of impedance are those computed by using equation (2) . 


DATA REDUCTION 

Perturbed-Flow-Rate and Total-Pressure Calculation 

The outlet perturbed flow rate is computed by using the following equation, which is 
a rearranged form of equation (9) of reference 4: 


Q. 


out 


p 9 - P10 
R or + 5 2rfL 9-10 


(3) 


where Pg - Pjq is the perturbed pressure difference measured across the orifice plate, 
R Qr is the resistance of the orifice plate, and L is the inertance of the pipe between 
the two measuring stations. The resistance R Qr is determined from the steady -state 
orifice calibration curve. Perturbed total pressures are computed from 

p = p + -£-Q q (4) 

A 2 


where A is the cross-sectional area of the pipe at the pressure measuring station, and 
Q av is the steady -state flow rate. In equations (3) and (4), the basic measured quanti- 
ties are the perturbed static pressures p. 


Impedance Calculation 

The pump impedance was computed from the perturbed total pressures and the per- 
turbed outlet flow rate by using equation (2) . This method is correct only for noncavita- 
ting data. The impedance is expected to be a complex number, with the imaginary part 
due to inertance. Thus, impedance can be written in the form 

P R- p s 

Z 6-8 = — q = R 6-8 + ^ X 6-8 = R 6-8 + ^ 2 , 7 fIj 6-8 

*^out 

where Z denotes impedance, R resistance, X reactance, and L inertance. The 
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symbols P and Q ref er to perturbed pressure and flow rate . Total pressures are 
used in the calculation as it is shown in appendix B that this is the correct choice for the 
general case of nonzero mean flow. The subscripts 6 and 8 refer to the measure- 
ment locations as shown in figure 1(b), and Zg_g refers to the impedance between these 
two stations. Zg_g is the combined impedance of the impeller and the volute. Two 
other impedances, Zg_g and Z g_g are also used; they are defined similarly . Zg_g 
is the impedance of the conical diffuser, and Zg_g is the impedance of the whole pump. 

Gain Calculation 

In the method of reference 4, pump gain is computed from two sets of data (the 
open-orifice data and the blocked -orifice data) taken at the same steady-state conditions 
but with different exit impedances in the test system. For the open orifice, equation (1) 
becomes 


P D,0 = (m + 1)P S,0 " ZQ 0 (6) 

(The subscript out has been omitted from Q ou f There is no chance of confusion, since 
only outlet flows appear in the equations.) For the blocked orifice, 

P D,B = ( m + 1 ^ P S,B " Z< ^B 

The pump gain m is assumed to be the same in both equations, so is the pump impe- 
dance Z. Solving the two equations for gain, 

, p d,o q b" p d,bQo , oX 

p s,o q b _p s,b q o 

This equation gives gain entirely in terms of measured quantities . 


RESULTS AND DISCUSSION 

All the data shown in this report were taken at the noncavitating inlet pressure of 
4. 83x1 0^ newtons per square meter gage (70 psig). (In ref. 5, the first effects of cavi- 

5 

tation on performance were observed at an inlet pressure of approximately 3. 79x10 
«> 

N/m gage (55 psig).) In order to study the effect of frequency on the pump dynamic 
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performance, data were taken at selected frequencies ranging from 2 to 52 hertz. The 
effect of flow rate was also investigated; the dynamic performance data were taken at 
three steady-state flow coefficients (0.442, 0.497, and 0.536) in the high-flow part of the 
pump operating range. Changes in flow rate result in changes in the impeller steady - 
state flow distribution, as discussed in reference 5. At higher flows, the flow rate 
through the impeller passage in front of the volute tongue is considerably higher than 
through the passages behind the tongue. Near the design flow coefficient of 0.363, there 
is approximately equal flow through all the passages. Some of the trends in the dynamic 
response of the pump can be explained to be the result of the changes in the steady-state 
flow distributions . 


Steady-State Performance 

The steady -state performance of the pump is shown in figure 5. Figure 5(a) shows 
the combined performance of the impeller and the volute, as measured just upstream of 
the conical diffuser. Figure 5(b) shows the performance of the complete pump (meas- 
ured at the diffuser outlet), and figure 5(c) the performance of the conical diffuser. It 
would be more conventional to show the diffuser performance in terms of a loss coeffi- 
cient; however, the head rise coefficient is used so that figure 5(c) conforms to the 
other two figures. The steady -state curves are of interest in dynamic studies since 
theoretically as the frequency approaches zero, the resistances should approach the 
slopes of the curves expressed in appropriate units (N-sec/m^, or lbf-sec/ft^). The 
steady-state pressure measurements were made at the same measurement stations 
(6, 8, and 9 of fig. 1(b)) as the dynamic measurements. Total pressure was measured 
directly at measurement station 6, but only static pressures were measured at stations 8 
and 9. The total pressures at the latter stations were obtained by adding the velocity 
head, using a velocity computed by dividing the overall flow rate by the flow area at the 
respective measurement station. The measured pressures at stations 6 and 8 agreed 
closely with those of reference 5 (in ref. 5, measurement station 8 is designated 5c). 

This comparison can not be made for station 9, which was not included in the reference 5 
measurements. 

There are some discrepancies in the open-orifice data. Four of the open-orifice 
data points of figure 5(b), taken on one day, agree with the bio eked -orifice data. Three 
of the open-orifice data points, all taken on another particular day, are all low. Unfor- 
tunately, the pressure at measurement station 8 was not recorded for the four good data 
points, so only the three questionable data points appear in figures 5(a) and (c). 
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Perturbed Flow Rates and Total Pressures 


The magnitudes of the perturbed pressures and flow rates are shown in figure 6 for 
a typical run. The greatest perturbed flow rate was about 4 percent of the steady-state 
flow rate (0. 0388 m 3 /see, or 1.37 ft 3 /sec) . The peaks at 17 hertz are due to changes 
in the acoustic impedance of the inlet line. Above 30 hertz the data are greatly affected 
by the longitudinal vibration of the inlet line, which had a resonant frequency of 34 hertz. 
This vibration is the cause of the differences between the flow out of the inlet line Qg 
and the flow into the pump Q^ n - The result is a low magnitude of Q- n which, in turn, 
affects all the pump perturbed pressures and flow rates . The impedance data (figs . 7 
to 10) are computed from ratios of the measured pressures and flows. Thus, the ups 
and downs of figure 6 do not appear in the impedance data. There is, however, an indi- 
rect effect. The measurement accuracy is poor when the measured quantities get too 
small. This explains the increase in the scatter of the impedance data which occurs in 
the frequency range of 30 to 40 hertz. 

Some analysis of the perturbed pressure and flow data was done in reference 4. 
Several checks were made on the consistancy of the perturbed flow measurements. The 
amount of data scatter in the compliance flow suggested an accuracy of no worse than 
2 percent in the inlet and outlet flow values. The noncavitating pump compliance was 
found to be large enough to cause significant differences to occur between the inlet and 
outlet perturbed flows. (A compliance value of 2.66x10”^^ m^/N (4. 5xl0 -3 ft^/lbf) was 
determined from the dynamic data.) Compliance is not a major concern in the present 
report, since all the results shown are based on the outlet perturbed flow. However, a 
complete representation of the pump dynamic response would have to include a compli- 
ance term to give the relation between inlet perturbed flow and outlet perturbed flow. 


Pump Impedance 

Three different impedances were computed; they are identified by their subscripts. 
Zg_g denotes the combined impedance of the impeller and the volute, Zg g the impe- 
dance of the diffuser, and Zg g the impedance of the whole pump. The impedance val- 
ues were computed by using equation (2) . Since this equation was derived assuming zero 
gain, it is valid only for noncavitating flow 

Impedance data . - The general trends of impedance with frequency are depicted in 
figure 7. The impedances Zg_g and Zg_g are shown for the open-orifice data at two 
flow coefficients (0.442 and 0. 536). (The blocked-orifice data have the same general 
trends as the open-orifice data.) The impedances are presented in terms of their real 
parts (resistance R) and their imaginary parts (reactance X) . All the reactances shown 
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are approximately directly proportional to frequency. This shows that they are caused 
by inertance, as anticipated. (Inertance is a parameter in hydraulic systems analogous 
to electrical inductance . ) A significant point is the relatively low frequency at which 
reactance becomes important. In fact, the reactance Xg_g becomes as large as the re 
sistance Rg_g at approximately 9 hertz at the high flow and 3 hertz at the low flow. 

In figures 8 to 10, the data are shown in terms of resistance and inertance. While 
figure 7 gives a better picture of the trends of pump response with frequency, figures 8 
to 10 are easier to analyze. This is because inertance is a convenient parameter to use 
in the electrical circuit analogy of a hydraulic system. All the noncavitating dynamic 
performance data are presented in figures 8 to 10, both the open-orifice and blocked- 
orifice data at all three flow coefficients. The dynamic representation of the whole 
pump, which is the item of interest in pogo studies, is given by Rg_ g and Lg_ g of 
figures 8(a) and (b), respectively. But in explaining the results, it is necessary to look 
at the separate contributions of the impeller -volute combination (Rg g and Lg g of figs 
9(a) and (b), respectively) and of the conical diffuser (Rg_g and Lg_g of figs. 10(a; 
and (b)). 

Data consistency . - One indication of data consistency is how well the impedance 
data of figures 8 to 10 agree for the two different outlet orifice configurations . The two 
configurations are identified as the blocked orifice and the open orifice . Theoretically , 
the data of figures 8 to 10 should not be affected by a change in outlet impedance, be- 
cause gain is assumed to be zero. The differences appearing in the figures should 
therefore be due to measurement errors and errors in setting the operating point (the 
same steady-state conditions could not be duplicated exactly for two different runs). 

In all the figures, there is good agreement between the open-orifice data and the 
blocked -orifice data. For most of the figures, no differences other than data scatter 
can be seen. However, for the inertance plots (figs. 8(b), 9(b), and 10(b)) at a pump 
flow coefficient Tp of 0. 497, there is a systematic difference between the two types of 
data of about 7 percent. The same percentage difference occurs in all three plots. 

This result, if it is due to measurement errors, would be caused by an error in Q out 
rather than by an error in one of the pressures p g , p g , or p g . An error in one of the 
pressures would affect only two of the three plots . The explanation of an error in Q Qut 
could be either an error in the measurement of the perturbed pressure across the outlet 
orifice plate, or an error in the steady-state flow measurement. (The latter would give 
an incorrect value of the resistance R Qr appearing in eq. (3) .) 

A few tests were run with a vibration aosorber added to the system to reduce the 
vibration of the inlet line. The results of these tests are shown as tailed data points in 
figures 8 to 10. There is no significant difference from the rest of the data. This indi- 
cates, as discussed in reference 4, that the effect of inlet line vibrations was correctly 
accounted for in the method used for reducing the data. 
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Impedance is a complex number , and its computed value will have an error in mag- 
nitude and an error in phase angle. At high frequencies, the reactive part of the impe- 
dance is much greater than the resistive part; consequently, an error in phase angle 
gives a comparatively large error in resistance. At low frequencies the opposite is 
true; so an error in phase angle gives a large error in inertance. In these cases, the 
individual resistance or inertance plot may show an increase in data scatter which does 
not indicate poorer measurement accuracy . 

Steady-state resistance values . - Theoretically, the resistance values at zero fre- 
quency can be calculated from the slopes of the steady-state head-flow curves. When 
dynamic data are not available, these steady -state resistance values can be used in an 
approximate analysis of the transient response of a pump. The steady -state resistance 
values were computed from the curves of figure 5. Because other curves could be 
drawn which appear to fit the data points equally well but have considerably different 
slopes, these resistance values are not expected to be very accurate. The slope of a 
performance curve will not be well defined unless the data points are closely spaced and 
unusually accurate. The steady -state resistance values are shown on figures 8(a), 9(a), 
and 10(a) by the symbol Rgg. The agreement with the dynamic resistance values ex- 
trapolated to zero frequency is probably within the accuracy to which Rgg can be deter- 
mined . 

Resistance and inertance trends with frequency . - Figure 9 shows that as frequency 
is increased, Rg_g increases while Lg_g decreases. The greatest changes occur at 
low frequencies; above roughly 20 hertz, changes are small. 

These trends can be explained by considering the effect of inertance on the pattern 
of perturbed flow through the pump passages . The pump internal geometry is shown in 
figure 11. At very low frequencies, the division of flow among the impeller passages is 
nearly the same as the steady -state flow division. This division tends to minimize the 
overall resistance. At high frequencies, inertance becomes important. The flow from 
the passages in front of the diffuser tongue takes a short, low-inertance path through the 
volute, as shown by the streamline labeled A in figure 11 . The flow from passages be- 
hind the tongue take long, high -inertance paths, as shown by streamline B. The high- 
frequency, perturbed flow travels more easily along the short, low-inertance paths. 
Thus, most of the flow will choose these paths, resulting in a flow distribution which 
tends to minimize the overall inertance, but which increases the overall resistance. 

This reasoning indicates that the resistance would tend to increase and the inertance to 
decrease with increasing frequency . 

At a high enough frequency, the flow distribution is determined mainly by inertance, 
with resistance having a relatively small effect. In such a case, a further increase in 
frequency will not change the flow distribution very much . This would explain why the 
changes in the data are small above 20 hertz. This high-frequency regime might be 
called inertance controlled. 
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A change in the steady-state flow rate appears to affect the inertance values L g_g» 
Lg_ 9 , and L g _ g more at low frequencies than at high frequencies. This also can be 
explained as a result of the flow distribution inside the pump. Changing the steady- 
state flow rate changes the resistances of the impeller passages. Whereas a change in 
frequency affects the perturbed flow distribution through changes in the inertive reac- 
tances, a change in flow rate affects the flow distribution through changes in resistance. 
It is the relative sizes of the various resistances and reactances in the impeller and 
volute which determine the flow distribution. As mentioned before, the effects of iner- 
tance overshadow the effects of resistance at high frequencies. Consequently, the 
changes in the impeller resistances affect the flow distribution more at low frequencies 
than at high frequencies. One result is that the overall inertance values, which are 
affected by the flow distribution, would tend to show greater differences at low frequen- 
cies than at high frequencies . 

At high frequencies, the differences become less than the data scatter, and the 
inertances approach common values which are unaffected by flow rate. In evaluating 
the inertance data at the flow coefficient of 0.497, it is useful to compare the high- 
frequency inertances to these common inertance values. At this flow coefficient, a sys- 
tematic difference between the open -orifice and blocked -orifice inertances was noted in 
the section Data consistency . The present comparison shows that this difference is 
probably due to errors in the blocked-orifice data, since the open-orifice data agree 
better with the common inertance values. 

At a flow coefficient of 0. 442, the diffuser inertance Lg_ g is constant within the 
accuracy of the data (fig. 10(b)). This shows that the diffuser is acting as a simple in- 
ertance element. The inertance can theoretically be calculated from the passage geom- 
etry by using the formula 



(9) 


where d l denotes an increment of passage length, A the cross-sectional area at posi- 
tion l, and B and C the beginning and end of the passage. The inertance computed by 
using this formula is shown as a solid line in figure 10(b). The agreement with the data 
at a flow coefficient of 0.442 is very good. 

At the other flow coefficients, Lg g deviates from the simple inertance model at 
low frequencies. It does not seem possible for the inertance to actually decrease below 
the theoretical value, as indicated. The average inertance would increase for any 
change in the flow distribution from the uniform flow assumed in deriving equation (9) . 
Therefore, it is suspected that the pressure measured at the diffuser inlet is not repre- 


12 



sentative of the average pressure at that location. At high flow coefficients, the flow 
entering the diffuser has high velocities near the tongue and low velocities near the outer 
casing. The static -pressure tap at measuring station 8 is located near the outer radius; 
so for nonuniform flow distributions , the measured pressure might differ significantly 
from the average pressure . At measuring station 9 (the diffuser outlet) , the velocity 
head is much smaller than at measuring station 8; furthermore, a nonuniform flow dis- 
tribution would probably smooth out by the time it passed through the diffuser. Conse- 
quently, the measurements at station 9 are regarded as more reliable than those at sta- 
tion 8, and the Rg_ g and Lg_ g values are considered to be more accurate than Rg_ 8 , 

R 8-9> L 6-8’ and L 8-9’ 

At low frequencies, Rg_ g is nearly constant for the flow coefficient of 0.442, but 
decreases with increasing frequency for the other flow coefficients. The diffuser resis- 
tance is basically caused by losses. One possible loss mechanism which would explain 
the observed trends is separated flow past the volute tongue. (Resistance is not propor- 
tional to loss, but to the derivative of loss with flow rate. But the resistance tends to be 
high when losses are high) . At high flows, the radial velocities near the tongue become 
much higher than their design values, so that the tongue operates far from its design in- 
cidence angle. This would cause the losses to increase with flow rate at low frequen- 
cies. Studies of airfoils in unsteady flow have shown that the tendency to stall decreases 
with frequency, so the losses would decrease with frequency. At high frequencies, the 
trends are different. Rg_ g tends to increase with frequency at all three flow coeffi- 
cients, and does not appear to be affected by flow rate. Since the perturbed flow distri- 
bution is determined mainly by inertance at high frequencies, the overall flow rate has 
only a minor effect on the data . 

To summarize, at the flow coefficient of 0. 442, which is the one closest to the de- 
sign flow coefficient, the dynamic behavior of the diffuser is simple and straightforward. 
The inertance is constant and can be computed from the diffuser geometry, while the re- 
sistance is constant (below 30 Hz) and equal to the steady -state value Rgg. However, at 
the two higher flow coefficients, the low-irequency trends are not so predictable. The 
reason for this may be the effects of the nonuniform flow distribution entering the dif- 
fuser . 


Pump Dynamic Performance Model 

A dynamic flow model can be developed for the pump by representing its component 
parts in terms of resistance and inertance. It has already been shown' that the conical 
diffuser follows the simple, constant -resistance, constant-inertance model near design 
flow. The pump impeller ana volute are a much more complicated case than the conical 
diffuser was, as there are 28 separate impeller passages, each emptying into the volute 
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at a different place. Each passage can be represented by an appropriate resistance and 
inertance in an electric circuit model of the pump. The impeller passages are treated 
as if they were not moving. It is shown in appendix B that for a radial-vaned impeller 
the inertance equation is the same whether the blades are moving or stationary. This is 
not true for a mixed -flow impeller. 

Electric circuit model of pump . - The pump passages are shown schematically in 
figure 11, and the electrical circuit model of the pump is shown in figure 12. Each of the 
28 impeller passages is represented by an inertance Lj and a resistance Rj. The vo- 
lute is represented by the circuit elements R Ai and L Ai , with i ranging from 1 to 27. 
The elements R^ and represent the resistance and inertance experienced by the 

volute flow in going from the outlet of one impeller passage to the outlet of the next pas- 
sage. The elements R A2 g and L A28 re P resent conical diffuser. Measurement 
station 9, located at the outlet of the diffuser, corresponds to the point having voltage 
E A29 * n c * rcu rt- Measurement station 8 is just before the diffuser inlet, the point 
having voltage E^ 2 g in the circuit. 

The representation of the flow passages by resistances and inertances is based on 
the assumption of a one -dimensional flow. The model requires that the flow conditions 
at any ©-coordinate of the volute be described by a single value of pressure and a single 
value of flow rate. Because of this requirement, the effects of radial pressure and flow 
variations can not be accounted for. The flow leaving each impeller passage is assumed 
to mix immediately with the flow already in the volute, whereas the actual case is more 
complicated. The model is intended to be only a relatively simple, approximate predic- 
tor of the pump dynamic response. 

The network of figure 12 was solved proceeding from left to right. Starting with a 
unit voltage for E A ^ , the current through Rj and Lj (representing the first impeller 
passage) is computed. This current passes through R A j and L A1 (representing the 
volute impedance), allowing the voltage difference between E^ and E A j to be de- 
termined. From E A2 , the current in the second branch of the circuit is determined. 
Continuing this calculation through all 28 branches of the circuit would be tedious, so the 
solution was programmed on a digital computer. The important results in this calcula- 
tion are the overall impedances. The impedance Zg_g, representing the response of the 
impeller and the volute, is equal to the voltage E A2 g divided by the total current (the 
current through R A2 g and L A2 g). The impedance Zg_g, representing the response of 
the entire pump, is equal to E A2 g divided by the total current. 

One effect of a change in frequency is to change the distribution of current in the net- 
work. The current distributions are compared in figure 13, which shows the magnitudes 
of the currents in the branches of the circuit which represent the impeller passages. 

The results are normalized by dividing the currents by the total current which flows 
through R^28 and L A28* The current decreases with increasing frequency for a 
constant overall voltage, because of the increase in the inertive reactances. At higher 
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frequencies, most of the current flows through the higher numbered branches of the net- 
work. In hydraulic terms, this means that the perturbed flow is concentrated in the pas- 
sages just before the tongue, as was mentioned previously. It is these changes in the 
perturbed flow distribution which are responsible for the changes in the overall resis- 
tance and inertance with frequency . 

The current distributions of figure 13 are for a case where the resistances from 
R A1 to R^27 were zero (the same case as was used to compute the zero-R^ curve of 
fig. 14). The zero-R^j case was chosen for this example since it best illustrates what 
actually happens in the pump as frequency changes. This case had equal flow in all pas- 
sages at zero frequency, which approximates the actual steady -state flow distribution 
close to design flow. Using nonzero R^ values resulted in unrealistic flow distribu- 
tions at zero frequency . 

Selection of resistance and inertance values . - To use the model requires a know- 
ledge of what the values of Rj, Lj, R^, and L Ai should be. While the inertance val- 
ues can be computed from the pump geometry, there was some question as to the best 
choice for the resistance values. Several choices of resistance values were tried, and 
the results were compared to the measured impedance data. What was desired was a 
relatively simple method of choosing the resistance values which would give a reasonably 
good data fit. How the inertance and resistance values were actually determined is dis- 
cussed in detail in the following paragraphs . Once a method is developed for determining 
all the parameters in the model, it can hopefully be applied to other pumps for which dy- 
namic data are not available. 

The inertance values were calculated from flow passage geometry by using equa- 
tion (9). The diffuser inertance 1^23 was 152 900 newton-second^ per meter (90.4 
lbf-sec^/ft^); this is the theoretical inertance value shown in figure 10(b). The impeller 
passage inertance was 720 000 newton-second per meter 3 (426 lbf-sec /ft 3 ) , while 
the volute inertances ranged from an value of 222 000 (131) to the value of 

17 500 (10. 36) . Even though the impeller passages were of different lengths due to the 
use of splitter blades, a single inertance value, based on the length of the main blade, 
was used for all. The inertance Lj also includes the inertance of a vaneless diffuser 
section, a narrow (0.635 cm by 2.54 cm, or 1/4 in. by 1 in.) passage above which the 
passage widens out to form the volute. The model requires that the flow be regarded as 
independent of the volute flow to a certain point, and completely mixed with the volute 
flow beyond that point. The flow out of a pump passage will actually take a little time to 
conform to the volute flow; so the dividing point was chosen to be somewhat beyond the 
impeller outlet. The L Ai values were calculated from the geometry of the volute (not 
including the vaneless diffuser section). The inertance values are inversely propor- 
tional to cross-sectional area, and are consequently high in the narrow part of the volute 
just behind the tongue, and low just in front of the tongue. 
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The values of Rj and R Ai can be chosen to match a given overall resistance value 
at a specified frequency. However, there are always many combinations of R^ and 
R Ai values which will do that. The effects of changing the relative sizes of the R Ai 
values compared to R^ are shown in figure 14. An estimate was made of the resis- 
tances occurring in the impeller and the volute, with the narrowest part of the volute 


assumed to have the highest R Ai 


"first estimate of R Ai 


values. " 


values. This case is identified in figure 14 as the 
For the second case (reduced R Ai values) the ratios 


R Ai /Ri were one -half those of the first case. For the third case (zero R Ai values), 


all R Ai values except R A2 g diffuser resistance) are zero. The Rj and R Ai 
values were scaled to give the same overall resistance Rg_g of 6. 26X10 6 N-sec/m 5 
(3700 lbf-sec/ft^) at 25 hertz. This resistance was chosen to match the data for the flow 


coefficient of 0. 442. 


The effects of changing the ratios R^i^l were important only at low frequencies. 
The zero-R A ^ case had the greatest changes in both Rg_g and Lg_g between low fre- 
quencies and high frequencies. This case also fit the experimental data best at the flow 


coefficient (0.442) closest to the design flow, although it still underestimated the differ- 


ence between the low-frequency and high-frequency Rg_g values. Because it is also 
the simplest case, the zero-R Ai case was selected as the best choice to use with the 
model . 


However , there is one point which should be considered since it casts a little doubt 
on the validity of the results shown in figure 14. The cases computed with nonzero R A ^ 
values had highly unlikely flow distributions. In one instance, the zero-frequency flow 
rate through passage 28 was over 100 times as great as the flow rate through passage 1. 
The question is whether these unrealistic flow splits jeopardize the results for the over- 
all R and L values. 


In order to see how a more reasonable flow distribution would affect the results, a 
second model was tried. The significant difference between the two models was that the 
one -dimensional -flow assumption of the first model, which made it impossible to con- 
sider the effects of radial pressure gradients in the volute, was discarded. In the sec- 
ond model, the flow from each impeller passage was represented by a separate stream- 
line in the volute; this allowed the radial pressure variations to be considered. The flow 
rates in the second model were specified to be equal in all impeller passages at zero 
frequency. Although this is a much more realistic flow distribution than that of the first 
model, the overall R and L values computed from the two models were not very dif- 
ferent. Because the items of interest are the overall R and L values and not the indi- 


vidual flow rates, the second model was discarded in iavor of the simpler first model. 

Comparison of model results to impedance data . - The results using the model are 
shown as solid lines in figures 8 and 9. The R Ai values (except R A2 g) were assumed 
to be zero, and Rj was chosen to fit the Rg_g data. To test validity of the model, Rj 
was chosen so that the overall resistance fit the data in the middle of the frequency 
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range, 25 hertz. The diffuser resistance R A2 g was chosen to fit the Rg_ g data. Dif- 
ferent values of Rj and R A2 g were required at each different flow coefficient. The 
model gave reasonably good curve fits for the Rg_g and Lg_g data of figure 9. It also 
agreed with the data for the entire pump (fig. 8) at the flow coefficient of 0. 442. How- 
ever, it did not agree with the low-frequency results at the other two flows. The lack of 
agreement is due to the deviation of the diffuser parameters Rg_g and Lg_g from the 
values assumed in the model. In general, the low-frequency results for these param- 
eters are hard to explain and may be in error, as previously discussed. More study is 
needed on how Lg_g and Rg_g vary with frequency. 

Application of model using steady -state data . - The model would not be very useful 
unless it could be used when dynamic data are not available. In this case, a different 
method of determining Rj and R A2 g * s necessary. The resistance Rj can be ob- 
tained from Rgg, which is proportional to the slope of the steady -state head -flow curve. 
For the zero-R Ai case, Rj is equal to Rg g times the number of impeller passages, 
since for this case the circuit of figure 12 becomes a simple parallel circuit at zero fre- 
quency. The resistance R A2 g is simply set equal to its steady-state value. Since the 
inertance values are computed from pump geometry, this means that the dynamic per- 
formance of a pump can be computed when only steady-state data are available. 

The results based on the steady -state resistance values are shown as dashed lines in 
figures 8 and 9. The agreement with the Rg_g and Rg_g data is best at the flow coef- 
ficient (0.442) nearest the design flow. An inaccuracy in the steady -state Rg_g value 
might account for most of the differences between the model results and the data. At the 
higher flow coefficients, the model overestimates the resistance values. In figure 9(a); 
an inaccuracy either in the steady -state R fi _g value or in the model at zero frequency 
would throw the results off at all frequencies. At the two higher flow coefficients, the 
inaccuracies in the model (as judged from the solid lines) are at least as important as the 
inaccuracies in Rgg- The inaccuracies in the Rg_g values used in the model are just 
as great as the inaccuracies in Rg_g at the higher flow coefficients; both affect the 
Rg _9 comparison of figure 9(a). The assumption that Rg_g is constant is a poor one 
at the higher flows . 

The agreement with the inertance data is better than with the resistance data, al- 
though the low-frequency Lg_g agreement is affected by the discrepancies in Lg_ g . 

The agreement of the model results with both the Lg_g and the Lg_g data is best at 
high frequencies. Furthermore, the inertive component of impedance becomes large 
compared to the resistance at high frequencies . This means that the model predicts the 
magnitude and phase of Zg_g and Z g _g more accurately at high frequencies. 

In summary, the flow model developed herein gave a reasonable prediction of the 
trends of pump impedance with frequency. It gave the best results at high frequencies 
and at the flow coefficient of 0. 442. The results are the most predictable at this flow 
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coefficient presumably because it is the closest to design flow, and consequently comes 
closest to having a uniform flow distribution in the impeller and volute. Obviously, be- 
fore the general usefulness of the model can be assessed, it must be applied to other 
pump configurations. However, the results for this one set of pump data have been en- 
couraging, particularly since the model provided reasonable estimates of pump impe- 
dance over a wide range of frequencies based only on pump geometry and steady -state 
flow measurements. 


Pump Gain 

Gain is generally assumed to be zero under noncavitating conditions. Since only a 
minute amount of tip vortex cavitation was present in the tests, the gain is expected to 
be very small. On the other hand, the method developed in reference 4 to determine 
gain (eq. (8)) is very sensitive to measurement errors. Thus the calculation of gain can 
be regarded as a check of the accuracy of the method of reference 4 as well as a check 
of the zero-gain assumption. The significance of the results can be judged with the help 
of an error analysis. The error analysis shows how much inaccuracy to expect in the 
computed gain values as a consequence of normal measurement errors. 

The values of gain computed for the entire pump are shown in figure 15 for all three 
flow coefficients. The discharge pressure P D corresponds to Pg, the suction pres- 
sure Pg to Pg. (Gain can also be computed for the impeller -volute combination, using 
Pg for Pp. The Pg results are shown here because the measurement of Pg is con- 
sidered more reliable than the measurement of Pg.) The gain values are seen to be 


TABLE I. - PROBABLE ERROR 
IN GAIN 


Frequency, 

Hz 

Root-mean-square 
error in gain for 
flow coefficient of - 

0.442 

0. 536 

2 

0.13 

0.41 

8 

.20 

.19 

14 

.28 

.26 

20 

.38 

.40 

30 

.55 

.39 

36 

1.33 

1.35 

40 

1.29 

1.44 

45 

2.20 

1.66 
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close to zero at low frequencies, but to diverge at high frequencies. 

In order to see whether the results are consistent with the assumption that gain is 
zero, an error analysis was done. The details of the error analysis are given in appen- 
dix C. The results are shown in table I. 

In general, the data points in figure 15 differ from zero by roughly the estimated 
error. The gain gets large just where the error analysis says it should, above 30 hertz. 
These results provide evidence that the observed trends in gain are due primarily to 
measurement errors, and that the error analysis provides an indication of the accuracy 
of the gain calculation method of reference 4. 

Gain can be checked in another way. The impedances calculated from the blocked- 
orifice and open-orifice data would not be expected to show close agreement if gain were 
significantly different from zero. Therefore, the effect gain would have on the value of 
Z nc will be examined. Equation (2) defines Z nc , the impedance computed assuming 
zero gain, as 


Z 


nc 



( 2 ) 


The suction pressure is given by 


Ps=-7T< P D + ZC W < 10 > 

m + 1 

This is simply a rearranged form of equation (1). The ratio of P D to Q out is set by 
the exit impedance Z 0 , 


P D " Z e Q out 

Substituting equations (10) and (11) into equation (2), 

Z nc =-i- ( Z- m Z e ) 
m + 1 

The difference between Z nc for the open orifice and the blocked orifice is given by 

Z nc,0 ' Z nc,B = ^7^( Z e,B- Z e,Q ) 


( 11 ) 


(12) 


(13) 
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Of course, measurement errors also produce differences between Z ~ and 
’ r nc , O 

Z t,. But if it is assumed that the data scatter is caused entirely by nonzero gain val- 

I1C/ j 1J 

ues, an upper limit to the magnitude of the gain can be computed. This was done for the 
flow coefficient of 0. 442 assuming a 5-percent difference, about the average data scat- 
ter, between Z nc q and Z nc g. (The real part of Z nc is Rg_ g , the imaginary part 
is 27 rfLg_ g ; the data scatter can be judged from fig. 8.) The exit impedance had a real 
part of about 6. 59x10^ N-sec/m^ (3900 lbf-sec/ft^) for the open-orifice case and about 
10. 65x10° (6300) for the blocked-orifice case. The imaginary part was the same for 
both cases. 

The magnitude of the gain which would explain the data scatter is shown in table n. 


TABLE II. - MAXIMUM GAIN CONSISTENT 
WITH THE AGREEMENTOF BLOCKED- 
ORIFICE AND OPEN-ORIFICE DATA 


Frequency, 

Hz 

Magnitude of maximum 
probable gain 

2 

0.09 

8 

.20 

20 

CM 

32 

.66 

45 

.92 


It can be said that, based on the data alone, the gain is not greater than these values. 
Furthermore, there is nothing in the data to suggest that the gain is other than the theo- 
retically assumed value of zero. 


SUMMARY OF RESULTS 

The dynamic performance of a radial -bladed centrifugal pump was measured using 
noncavitating inlet pressures. The results were analyzed in terms of resistance and 
inertance; this is normal procedure for a hydraulic system. Some of the results were 

1 . The reactance is important at all frequencies and is greater than the resistance 
except at relatively low frequencies (below 3 to 9 Hz, depending on the flow rate). 

2. The resistance values calculated from the steady-state performance curves 
agreed reasonably well with the dynamic resistance data extrapolated to zero frequency. 

3. Both resistance and inertance are frequency dependent. They change most rap- 
idly at low frequencies; above 20 hertz, the changes are small. For the impeller -volute 
combination, the resistance increases with frequency, while the inertance decreases. 
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4. At the flow rate nearest the design flow, the inertance data for the outlet diffuser 
agreed very well with the inertance computed from the diffuser geometry. The diffuser 
resistance was constant (below 30 Hz) and equal to the steady -state value. At higher 
flows, some unexpected trends occurred at low frequencies, possibly caused by the non- 
uniform flow distribution at the diffuser inlet. 

A model was developed in order to explain the trends of resistance and inertance 
with frequency. In this model, the pump is represented by an electric circuit. The 
model does not require any dynamic data; inertances are computed from pump passage 
geometry, resistances from the steady-state performance curves. The results using the 
model showed good agreement with the data at the flow closest to the design flow. 

The agreement is not as good at the higher flows, presumably because the perform- 
ance is affected by the uneven flow distribution in the impeller and the volute. At these 
higher flows , the model gave reasonably accurate predictions of overall pump impedance 
at high frequencies, but not at low frequencies. The model overestimated the overall 
resistance values at almost all frequencies. This happened because the model assumed 
the individual passage resistances are constant, whereas they are actually frequency de- 
pendent. To improve the model results would require a better understanding of the 
losses due to off -design flow patterns. The inaccuracies in resistance are not critical at 
high frequencies, where the inertive reactance is the main component of pump impe- 
dance . The model results agreed with the inertance data better at high frequencies than 
at low frequencies. The discrepancies in the diffuser inertance at low frequencies and 
high flows also merit further investigation. 

The data were also used to calculate gain. This was done both as a check on the 
method developed for determining gain and as a check on the assumption of zero gain. 

The gain was found to be zero within the accuracy of the measurements as shown by an 
error analysis. The error analysis also showed that the gain calculation is very sensi- 
tive to measurement errors. 

Lewis Research Center , 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 23, 1971, 

128-31. 
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APPENDIX A 


SYMBOLS 

2 2 

A cross-sectional area, m (ft ) 

E voltage 

F blade force vector 

f perturbation frequency, Hz 

j imaginary unit, V^I 

L inertance, N-sec^/m^ (lbf-sec^/ft^) 

l length coordinate 

l direction of blade passage 

m gain (see eq. (1)) 

P perturbed total pressure, N/m^ (lbf/ft^) 

p perturbed static pressure, N/m (lbf/ft ) 

3 3 

Q perturbed flow rate, m /sec (ft /sec) 

•j 3 

CD steady -state flow rate, m /sec (ft /sec) 

clV 

R resistance, N-sec/m^ (lbf-sec/ft^) 

r radius vector 

t time, sec 

U wheel speed, m/sec (ft/sec) 

V velocity, m/sec (ft/sec) 

W relative velocity, m/sec (ft/sec) 

X reactance, N-sec/m^ (lbf -sec/ft ^) 

x length coordinate, m (ft) 

Z impedance, N-sec/m^ (lbf -sec/ft ^) 

A finite difference or, in appendix C, magnitude of a random error 

p density, kg/m^ (slug/ft’*) 

( p flow coefficient (ratio of average inlet flow velocity to blade tip speed) 

ip head -rise coefficient 

w rotor rotational speed, rad/sec 
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u>' 2ir times perturbation frequency 

Subscripts: 

A volute circuit component in fig. 12 

av average 

B blocked orifice 

D discharge (either measurement station 8 or station 9) 

e exit 

i branch of network, ranges from 1 to 28 

in inlet 

nc noncavitating 

0 open orifice 

or orifice 

out outlet (measurement station 9) 

R relative to moving coordinate system 

5 suction 

SS computed from steady -state data 

6 peripheral coordinate 

1 impeller passage circuit component in fig. 12 

1,2 beginning and end of flow passage, appendix B 

5,6, 8, measurement station locations, fig. 1 (A double subscript separated by a 
9,10 hyphen (as in Zg_g) refers to the flow passage between the respective mea- 

surement stations.) 

Superscript: 

perturbation quantity 
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APPENDIX B 


DERIVATION OF INERTANCE EQUATION 

The inertance equation is usually derived for a quiescent fluid with stationary 
boundaries . The first part of this appendix derives the inertance equation for the case 
of a flowing fluid. The second part considers the case where both the fluid and the 
boundaries are moving. 


Derivation of Inertance Equation for Stationary Boundaries 

The analogy between fluid inertance and electrical inductance is commonly used in 
the field of acoustics. The inertance equation, on which this analogy is based, is usually 
derived for the case of a disturbance acting on a fluid at rest. The following derivation 
considers the case of a disturbance superimposed on a fluid in motion. The result is an 
equation having the same form as the usual inertance equation, but using total pressure 
instead of static pressure. 

The equation of motion for a fluid, neglecting viscous terms, is (ref. 6, eq. 3(9), 
neglecting all forces except pressure) 


DV 

p — = -VP 
Dt 


(Bl) 


Applying this equation to a simple one -dimensional flow as shown in figure 16, 


'9V + v 3V\ = _ap 


at 


9x 


9x 


(B2) 


Making the following substitution 


p + V 2 =i£ +pV 5V 

9x 9x \ 2 / 9x 9x 


(B3) 


gives the equation in terms of total pressure 


9V _ _9P 
P 9t 9x 


(B4) 
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Now assume that V is composed of a steady-state component V ay and a perturba- 
tion component V - ', as given by the following equation: 

v = V av + Ve ja,lt (B5) 

In this equation, the phase differences in the perturbation term at different x -coordinates 
are neglected. This is a good approximation if the length of the passage is much shorter 
than one -quarter wavelength. A similar equation is assumed for P, 

P = P ay + P'e^’ 1 (B6) 

Substituting equations (B5) and (B6) into equation (B4) and separating out the time depen- 
dent terms yield 


• ’ Ti-’ dP’ 
joj pV 

dx 


(B7) 


Integrating over the length of the passage (from x^ to Xg) , 


- y* 2 dP’ = J 2 j co'pV' dx = jo/ f 
X 1 X 1 «/x. 


AV*p — 
A 


(B8) 


The mass flow pAV' can be assumed to be constant if the passage is much shorter than 
one-quarter wavelength. Equation (B8) can then be written as 


P - p' - 
*1 2 


t . . » /* X 2 dx 

2 = AV p I — 


(B9) 


This is the standard form of the inertance equation. The total pressure P' is analogous 
to voltage, the volumetric flow AV* to current, and the inertance p J' dx/A to induc- 
tance. 
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Derivation of Inertance Equation for Rotor Passages 

One form of the equation of motion for relative flow in a turbomachine is (ref. 6, 
eq. 7(39)) 


at 


+ - VP 
p 


R 


= W x 


(curl W + 2to) 


+ F 


(BIO) 


where the relative total pressure P R is given by (neglecting the gravitational force 
term) 


P 


R 


= P + 




(Bll) 


Assume that there are a large number of blades and relatively narrow passages . 

The relative velocity W is constrained to flow in the direction of the blade passages. 
This direction is designated as l in figure 17. Since the vector W x(curl W + 2to^ is 

perpendicular to W, it has no component in the l direction. Neglecting viscous forces, 
F represents the force exerted by the blades on the fluid. This force is normal to the 
blade surfaces, so it also has no component in the l direction. Consequently, the com- 
ponent of equation (BIO) in the l direction is simply 


3 rW SP 5 

p 

at di 


(B12) 


The quantities W and Pj^ can be assumed to be composed of a steady component 
and a perturbation component, as was done in the preceding section for V and P: 


W = w av + WV W 1 


P R “ P R,av + P R e ^ t 


(B13) 

(B14) 


Note that to and to' are not the same; to is the rate of rotation of the turbomachine in 
radians per second, while to' is 2n times the perturbation frequency . Substituting 
equations (B13) and (B14) into equation (B12) and keeping only the time -dependent terms, 
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(B15) 


joj’pW’ = - 



Integrating over the passage length yields 

J r 1 2 

1 i”' AW 'p 

h 

However, to analyze a system containing both rotating and stationary flow passages 
would require that this equation be given in terms of absolute pressures and velocities. 
For velocities, the perturbed component is the same whether the flow is relative or ab- 
solute. Expanding the equation 



V=W + ^o>Xrj 

in terms of steady and perturbed components (it is assumed that To is constant) 

V' + V'e^'* = W a _ + W'e ja,,t + (w x r) 

cLv aV \ / 

Equating the two time -dependent terms, 


(B17) 


(B18) 


V’ = W' (B19) 

For pressures , start with the relation between total pressure and relative total 
pressure 



" P R + 


,(uw e + u 2 ) 


(B20) 
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where U denotes the wheel speed which is the magnitude of aJ x r and W 0 is the com- 
ponent of W in the u xr direction. 

In terms of steady-state and perturbed components, 


p + t - p + P , 

F av + Fe -^R,av + F R' 


k ei "’‘ + p[ u ( w e,av - + u2 ] < B21 > 


Equating the perturbed components , 


P’=P^ + P UW’ 


(B22) 


Substituting equations (B22) and (B19) into equation (B16), 


P i-P u l v 9,l- P 2 + P u 2 v e, 2 " j "’ Av 'p/ T- 


i 


dl 


(B23) 


Unless there is turning upstream of the pump inlet, ^ will usually be zero. 
Furthermore, for a radial -vaned centrifugal pump, ^ ^ = 0. Thus, for the test 

pump of this report, both terms are zero, and equation (B23) becomes the same as 
the equation for stationary passages (eq. (B9)). 
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APPENDIX C 


ERROR ANALYSIS OF GAIN CALCULATION 

The effect of random measurement errors on the accuracy of the gain values can be 
determined by using the methods of reference 7. An equation can be derived giving the 
error in gain as a function of the errors in the individual measurements . Let f denote 
a function of the measured variables x^, Xg, . . . , x n - If the errors are assumed to be 
random errors, the probable error in f is given by (ref. 7, eq. (3)) 



where A denotes the magnitude of the respective error. 

The function to be considered in this analysis is the gain m as given in equation (8) 


m + 1 


p d,o q b p d,b q o 
p s,o q b _p s,b q o 


(C2) 


d(m ♦ 1) — L [<DP d 0 - NP g c )dQ B - <DP d b - NP S 


«B< D dP D,0 - N dP S,0> - dP D,B - N dP S,B>] < C3 > 


where N denotes the numerator of equation (C2) , and D the denominator . 

The perturbed flow rate Q is calculated from several measurements, one of which 
is Pp. So before equation (Cl) can be applied, the errors in Q due to errors in P^ 
must be separated out and combined with the other terms containing dP^. The two 
parts of dQ shall be written as 


dQ 



dP 


D 


+ dQ other 


Applying equation (Cl), the error equation for gain is 


(C4) 
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DQ b - (DP d>b - np S jB ) 


3P 


0,0 



DCi o " (dp d,o " NP s,o ) 


3Q 


B 


3P 


D.B 


2 4 

( ap d,b) + I n 9q 2 I ( ap s,b) 


(C5) 


Now consider the outlet flow rate Q more closely. It is calculated from 


i(p 9 - P 10 + P 9 . 10 > 

Q = 


2 AP + 

Q 




(C6) 


where Pg is the static perturbed pressure measured in front of the multihole orifice 
plate, Pjq the pressure measured behind the orifice plate, and Pg.^Q a directly mea- 
sured pressure difference across the orifice plate. Taking the derivative of this equa- 
tion and dividing by Q gives 


dQ dp 9 - dp 10 + dp 9-10 
Q p 9 ■ PlO + p 9-10 



(C7) 


The outlet line reactance w L g_io is sma ^ compared to the multihole orifice resistance 
2 Ap/q, and has been neglected. Measuring station 9 is the discharge measuring sta- 


tion; so pg is the static -pressure component of the discharge total pressure P D 
errors in the velocity head P B - Pg are small compared to the errors in Pgj so Aj 
is approximately the same as AP^. Consequently, the partial 3Q/3 P D is given by 


The 
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Q 


(C8) 


9Q Q 

9P D P 9 " P 10 + p 9-10 2p 9-10 

(Since P 9 _ 10 is nearly equal to p g - p 1Q , the denominator equals 2p g _ 1Q accurately 
enough for these calculations.) The other terms in equation (C7) combine to form 
AQ ot her ’ which is Siven by 

Hotherf ( ap lo) 2 * ( ap 9-lo) 2 + 

Q 2 1 2p 9-10 1 

Using equations (C5), (C8), and (C9), the error in gain can be calculated from esti- 
mates of the errors in the individual measurements. Two types of error were con- 
sidered in estimating the amount of error in the perturbed pressures: (1) a 48-newton- 

o 

per-square-meter (1-lbf/ft ) error representing the background noise due to random 
pressure fluctuations, and so forth-, and (2) an error of 2 percent of the measured value. 
These two errors were combined by using the sum-of-the-squares rule for random er- 

9 

rors. The 48-newton-per -square-meter (1-lbf/ft ) and 2 -per cent values are judgements 
based on experience with this type of measurement. The outlet orifice resistance 
2 Ap/q was assumed to have a 2-percent error. 

The results of the error analysis are shown in table I (RESULTS AND DISCUSSION 
section) of this report. A simplified form of equation (C5) was used in which it was as- 
sumed that N = D (i.e. , that the gain is zero). The error in gain was found to be much 
greater at high frequencies than at low frequencies. The values in table I are highly ap- 
proximate, and are intended to show mainly the order of magnitude of the errors. Be- 
cause the two terms in the denominator of equation (C2) nearly cancel out, a relatively 
small change in any of the variables in the denominator could significantly affect the sen- 
sitivity of the gain to measurement errors. Also, the measurement error estimates are 
of a highly approximate nature. 
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(a) General view of test facility. 

Figure 1. - Pump Pertu roations Test Facility, showing measuring stations. For analytical purposes, stations 1 to 5 and 6A are assumed to be attached to 
ground; stations 6 to 10 are assumed to be attached to the pump and pipe. 





Measurement station 9, 
downstream of diffuser outlet 


Conical diffuser 



Measurement station 6 



(b) Test pump, showing measurement stations. 
Figure 1. - Concluded. 
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Figure 3. - Radial-biaded centrifugal pump impeller. 
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Figure 4. - Multihole orifice plate. 
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Steady-state noncavitating performance of centrifugal pump. 
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Figure 6 . - Magnitude of perturbed pressures and flow rates. Flow coefficient, 0.442; open orifice 
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(a) Flow coefficient, 0.442. 


Figure 7. - 
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(b) Flow coefficient, 0.536. 


impedance - open -orifice data. 
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(a -3) Flow coefficient, 0.442. 


la) Resistance, R 6 _ 9 . 

Resistance and inertance across the complete pump. 
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(a-1) Flow coefficient, 0.536. 
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(a -2) Flow coefficient, 0.497. 
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(a-3) Flow coefficient, 0.442. 

(a) Resistance, R 6 _ 8 . 

Figure 9. - Resistance and inertance across impeller and volute. 
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(b-1) Flow coefficient, 0.536. 
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(b-2) Flow coefficient, 0.497. 







□ Blocked -orifice data 
O Open -orifice data 
d Vibration absorber used 


O 

P DO 

O ° 



(a-1) Flow coefficient, 0.536. 


-°R 

^, R ss 


K ss 


cP n ° □ 


° cp cP B ° 
□ IP o-, P-, 


□ 


Qd Rh 


(a-2) Flow coefficient, 0.497. 
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(a -3) Flow coefficient, 0.44Z 
(a) Resistance, R g _ 9 . 

Figure 10. - Resistance and inertance across conical diffuser. 
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(b) Inertance, L^_g. 

Figure 14. - Computed resistance and inertance across impeller and volute for different relative sizes of volute resistance 
to impeller resistance. 
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Figure 15. - Pump gain. 
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Figure 16. - One-dimensional flow. 



Figure 17. - Pump passage geometry. 
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